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Abstract 

CL ' We study a nonlinear porous medium type equation involving the infinity Laplacian operator. We 

^^ ' first consider the problem posed on a bounded domain and prove existence of maximal nonnegative 

viscosity solutions. Uniqueness is obtained for strictly positive solutions with Lipschitz in time data. We 

i-Q also describe the asymptotic behaviour for the Dirichlet problem in the class of maximal solutions. We 

then discuss the Cauchy problem posed in the whole space. As in the standard porous medium equation 
(PME), solutions which start with compact support exhibit a free boundary propagating with finite 
speed, but such propagation takes place only in the direction of the spatial gradient. The description 
of the asymptotic behaviour of the Cauchy Problem shows that the asymptotic profile and the rates of 

£NJ ■ convergence and propagation agree for large times with a one-dimensional PME. 

> 

^ : 1 Introduction 

r — . In order to understand modes of nonlinear diffusion with preferential propagation in some directions, we 

C^> ' consider here the following variation of the porous medium equation, 

o: 

(1.1) Pt = Aoo(p m ), 

where m > 1 and Aoo denotes the 1- homogeneous version of the infinity-Laplacian, 



& ' Aqo u := \Du 



- 2 E 



X ' x ■ cc ■ n 



This equation is related to the standard porous medium equation (PME): p t — A(/9 m ) = 0. It represents a 
diffusion of porous medium type, but taking place only in the direction of the spatial gradient of p. Moreover, 
if a solution of (|1.1[) is radial (i.e., u = u(\x\,t)), it corresponds exactly to a solution of a 1-d PME with 
spatial variable r = \x\. 

Here we examine the general non-radial case with the intention of proving that there exists a natural 
evolution process that still looks like a 1-d directional PME, at least for large times. However, in this case 
the diffusion chooses the direction of the gradient, that varies with space and time. The first difficulty in 
devising a theory is to define a suitable concept of generalized solution for this nonlinear equation, which is 
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formally parabolic but can be degenerate or singular at some points: Here we use a well-known trick of the 
PME theory, the so-called density-to-pressure transformation, and introduce the variable 

(1.2) u:= _!2_pm-l. 

TO — 1 

Proceeding formally, it is easy to see that, if everything is smooth, u must solve 

(1.3) u t = kuA ao u+ \Du\ 2 

with k = 77i — 1 > 0. We propose to work out the existence and uniqueness theory for the latter equation, 
translating back the results to equation (|1.1|) via (|1.2|l . We study such problems both in a bounded domain 
and in the whole space, and we establish regularity, finite propagation, and asymptotic behaviour under 
convenient assumptions on the data. 

1.1 Results on a bounded domain 

We couple equation (|1 .3[) with initial and boundary data: 

(1.4) u{x,t) = g(x,t) on F, 

where T = d p Q is the parabolic boundary if our domain Q = £1 x [0,T], fl c M. d . As mentioned above, if 
the problem is posed in a ball with radially symmetric data it is observed that the solutions are equivalent 
to solutions for the 1-d porous medium equation, for which there is a weak theory that gives existence and 
uniqueness. Such a theory is not available for more general solutions of the fully- nonlinear equation (jl.3p . 
Therefore, we introduce the concept of viscosity solution. Since the equation is both singular and degenerate, 
we adapt the modified notion of viscosity solution introduced in |CVj and [BVj and the viscosity method to 
tackle degenerate problems. Our first main result concerns solutions with strictly positive data. 

Theorem 1.1. Let fi be an open and bounded domain in M. d . For each g £ C(T), g ^ c, there exists a 
viscosity solution u £ C(Q) of (|1.3|) satisfying u — g on T . Moreover, u Js c and the modulus of continuity 
of u can be estimated in terms of the modulus of continuity of g; in particular, if g is Lipschitz, then so is 
u. If g is Lipschitz on the lateral portion ofT, then the solution is unique. 

Our main interest lies however in treating nonnegative data that may vanish somewhere, in the spirit of 
porous medium equations with finite speed of propagation. This is not so easy. If we do not assume the 
data to be strictly positive, we can still obtain a maximal viscosity solution. 

Theorem 1.2. Given g £ C(T) with g Js 0, there exists a viscosity solution of (|1.3|) . u £ C{Q), satisfying 
u = g on T , and such that if v is another viscosity solution of (|1.3p with the same initial and boundary data, 
then ii^ u. Moreover, if g is Lipschitz continuous, then u is also Lipschitz continuous. 

Maximal viscosity solutions are our choice of good class of solutions to work with. 

Let us outline the organization of the proofs of this and related results. In Section [2] we define viscosity 
sub- and supersolutions for this problem and prove a useful lemma to simplify testing a solution. 

In Section [3] we consider strictly positive solutions. We prove two comparison results: a "weak" and a 
"strong" version. The strong version is valid for functions which are Lipschitz on the lateral boundary, and 
is in fact just a regular comparison result. 

In Section |4] we compute some explicit solutions. Of particular interest are the Barenblatt functions, which 
are similarity solutions, and the traveling waves. 

For existence, the idea is to regularize the problem as is done in |JK] for the infinity Laplacian evolution, 
and obtain estimates independent of the approximation process. In Section [5] we prove the existence of 



maximal viscosity nonncgative solutions and give some properties of these solutions, in particular we prove 
a uniqueness result for this problem in a nice domain (slightly more than star shaped). 

Section |5] is the longest and contains the more technical proof of the main result, Theorem ll.il 

To complete the study in bounded domains, Section [7] discusses the large time behaviour of solutions of the 
Dirichlet problem, which is described in the following theorem. 

Theorem 1.3. Let u be the maximal solution of (|1.3|) - (|1.4[) with g(x,t) = for x E dil. Then, 

(1.5) lim tu(x,t) = U(x) = U n {x) = [F n (x)] m -\ 

t— >oo 

where U is a Lipschitz continuous function, positive in fi which vanishes on dfl. 

As an interesting consequence of our study here, we obtain existence of solution of the elliptic eigenvalue 
problem 

(1.6) -A 00 G = \G p , p<\. 



As in |CV| and jBVj . we approximate the data with a strictly positive function and then take the limit. This 
allows us to obtain the asymptotic behaviour of solutions for the Dirichlet problem in Theorem 1 1.3 1 as saying 
that u(x,t) approaches the so-called Friendly Giant, U{x,t) = t Fq(x), whose existence is well-known in 
the PME theory. Our results are also consistent with the asymptotic profile obtained by Laurengot and 
Stinner |LS| for the infinity heat equation. 

Remark. We will show that the function Fq is in fact a constant multiple of the solution of the elliptic 
eigenvalue problem in fl (|T.3|) . 



1.2 Results on the whole space 

In Section|8]we discuss the Cauchy problem, posed in the whole space for nonnegative solutions. We assume 
that the initial function is nonnegative and bounded. We prove existence of a maximal solution. We can not 
obtain a complete comparison result in this case, but for compactly supported data the results for maximal 
solutions are valid, namely the propagation of the free boundary and its regularity in the form of the following 
theorem. 

Theorem 1.4. Let us assume that po : R d — > R is continuous, nonnegative, bounded and compactly sup- 
ported. Ifu(x,t) is the maximal solution of the Cauchy problem (j8.ll) . then there exists a constant R > 

which depends on m,n and the initial data po, such that as t — > oo, p(x,t) = [ m ~ 1 ul satisfies 

(1.7) t^ \p(x,t) - f3 R (x,t)\ ^ 

uniformly in x £ W l , where /3r is the Barenblatt function defined in (|4.3j) . Moreover, we have convergence 
of the supports 

(1.8) B Rl(t) (0) C {x : p{x,t) > 0} C B fl2 ( t) (0) 

where R i (t)/t 1 ^ m+r > -> R = R(m,n,p ). 

The Barenblatt function p R (r,t) is explained in formula (|4.3j) and R is its radius at t = 1. There is to our 
knowledge no simple formula to express the dependence of the asymptotic constant R on the data pq. Such 
a difficulty is relatively frequent in problems in nonlinear mechanics, see |KP Vj . 

The paper ends with an appendix on special solutions. 



1.3 Extension 

A natural development of the main idea of this work leads to consider similar models of propagation in a 
preferential direction. One option is to use interpolation of the regular Laplacian with the infinity Laplacian. 
We propose as the simplest example the family of equations 



(1.9) Pt =L E (p m ), where L e = e A + (1 - e) A M 

with < e < 1. Since the standard p-Laplacian operator, 1 < p < oo, is defined as 



A p u = div(\Du\ p ~ 2 Du) = \Du\ p - 2 ^u x%Xj j 



(1.10) A p u = div{\Du\P-'Du) = \Du\p-' > u XiXj I % + (p 



i,3 



\Du\ \Du\ 



if we put e = l/(p — 1) we can also write this proposed model as 

(1.11) Pt = e\D{p m )\ 2 -PAp{p m ). 

Many of the results of the present paper apply, at least partially, to the interpolated family. In particular, 
we note that, at least for radial solutions, (|1.9|) is like a (PME) in dimension 1 + e(d— 1). We will not pursue 
such an analysis in the present work. 

Notations 

We consider equations defined on some subdomain of the whole Euclidean space-time S = M. d x E. For 
a point Pq = (xo,t) £ S, we say that U is a parabolic neighbourhood of Pq if Pq £ Q C U, where Q is a 
cylinder centered at Pq, that is, Q = B t (xq) x (t — r, to] for some r, t > 0. We denote by C + (D) the space 
of nonnegative continuous functions from D to R, whereas C 2,1 (D) will denote those functions which are 
twice differentiable in x and once in t. Whenever Q is a cylinder of the form Q = X [0, T], with H, C K d 
open, we denote its parabolic neighborhood by T = d p Q = {dtt x [0, T]) U (CI x {t = 0}). 

For a function u £ C 2,1 (D), Du and D 2 u will denote the spatial gradient and the d x d matrix of second 
derivatives of u, respectively, 

Du= (u Xl ,...,u Xd ), D 2 u= {{u XiXi )).. =1 . 

Given a symmetric d x d matrix with real coefficients A, we define 

A(A) := max (Aw) ■ w, X(A) := min (Aw) ■ w, 

in other words, the largest and smallest eigenvalues of A, respectively. We will also need an approximation 
function, j3 c , which is a smooth real function satisfying /3 c (z) = \ z \ if \ z \ ^ c and /3 c (z) ^ c/2 everywhere. 

2 Viscosity solutions 

Following [CVj and |BVj . we define viscosity solutions for the modified equation (|1.3|) . Note that because the 
equation is singular at points where the gradient of the function vanishes, the usual definition of viscosity 
solution needs to be adapted at the singular points. We adapt the definition from |CGG| . see also jJKj . 

Definition 2.1. Given u £ C + (Q), we say u is a nonnegative viscosity subsolution of (|1.3[) in Q if 
and only if for every Pq £ Q and every function <fi £ C 2,1 (Q) which touches u from above at Pq, the following 
holds at the point Pq : 

(jJt^k^A^cjy+lD^ 2 ifDc^^O, 

4> t ^ k(f)A(D 2 (t>) ifD<t> = 0. 



As in I JK wc can weaken the second condition. 



Lemma 2.2. Let the function u satisfy the following: given Pq £ Q and </> € C 2 ' l (Q) such that u — <f> has 
an absolute maximum, at the point Pq and u(Pq) = 4>{Pq), it follows that at the point Pq 



b t ^0 ifDcf> = 0, D 2 (j) = Q. 



Then u is a viscosity subsolution of (JTT3J) . 

Proof. Step 1. Assume u is not a viscosity subsolution of (|1.3[) but satisfies the condition of the lemma. 
Then there exist Pq — (xo,to) *= Q anc ^ ^ C 2,1 (Q) such that u — <f> has an absolute maximum at Pq, 
u(P ) = <t>(P ), D(j)(P ) = 0, D 2 <j>{P ) ± 0, and at P 



<t> t > kc/)A(D 
Let us define 



2, 



3 l~ -.14 3 u „|2 



Wj(x,t,y,s) :=u(x,t) ~4>{y,s) - -\x-y\ - -\t - s\ 

and let (xj,tj,yj,Sj) be a point of maximum for Wj in Q x Q. It is easy to see that (xj,tj,yj,Sj) — > 
(xo,to,Xo,to) as j -> oo. 

Step 2. Let us check that for j large enough we can not have Xj = y,-. Assume iej = j/j and let us define a 
new function 

e j (y, s ):=-i\x j -y\*- J -\t j -s\ 2 

Then 4> — 9j has a local minimum at the point (j/j, s^) and hence, at this point <^ t = 6j t t and _D 2 Js L* 2 ^ = 0. 
From our assumption, for j large enough wc have 

j(tj - Sj) = 0j,t(yj, Sj) = <t>t{Vj,Sj) > k(j){y j ,s j )K{D 2 (j){y ] ,s ] )) ^ 0. 

Similarly, with 

Q( x ,t):=l\x- yj \ 4 + J -\t- Sj \ 2 

the function u — Q has a maximum at (xj,tj) and both DQ and D 2 Q vanish at (xj,tj), because Xj = ijj. 
Therefore, Q + C satisfies the condition of the lemma and we conclude that 

j(tj -Sj) =Cj,t(xj,tj) < 0, 

a contradiction. 

Step 3. We have concluded that Xj ^ j/j for large j. Let us now check that this also leads to a contradic- 
tion. By our assumption, and using the continuity of F H A(D 2 <j>(P)), there exists e > such that in a 
neighbourhood of (x ,tg) 

e< cj> t -kcj>k{D 2 <i>)- \Dcj)\ 2 . 

Therefore, for large j, using the fact that (f> — 9j has a minimum at (jjj, Sj), 
e< (4>t~kcj) A oo 0-|L>0| 2 )(y J ,s J ) 
0(yj, s j) ,r>2 t 



< j(t 3 - Sj ) - * ,;, ..'la ^gjfa.gj) • fo - Vj), (xj - Vi)) 



flxj-Vjf. 



On the other hand, we can still apply the condition of the lemma to Q, but now DQj ^ 0. Using this and 
the fact that u — Q has a maximum at (xj , tj ) we have 



j(tj - Sj ) - j'\ X j - yj \ b sc k K 3 \' {D\j{ X j,tj) ■ (xj - Vi ), (xj - Vj )). 

\ x j yj\ 

Since D 2 (j(xj,tj) = —D 2 9j(yj,Sj), this and the above equation lead to a contradiction for large j. □ 



Definition 2.3. Given u £ C + (Q), we say u is a nonnegative classical free- boundary solution of 
(|1.3[) if and only if: 

(i) on the positivity set V(u) = {P £ Q \ u(P) > 0}, the function u is smooth and solves (|1.3[) in the 
classical sense; 

(ii) the boundary of the positivity set T = dV(u) C\Q is a smooth hypersurface and u £ C 2 ' l (V(u) n T); 
(Hi) on the hypersurface T we have 

(2.1) u n = \Du\, 

where o~ n denotes the normal speed of boundary T. 
If instead of condition (i) above we have 

u t ^ ku Aqo u + \Du\ , 
in the positivity set and instead of (|2.1[) we have 

a n sC \Du\, 
then we say that u is a classical free-boundary subsolution. If in Definition ^. 31 wc impose the extra condition 
O) \Du\ ^ on T, 

then we say that u is a classical moving free-boundary solution. Replacing all the inequalities with ^ by 
inequalities with ^ we define classical (moving) free-boundary supersolution. 

Given two functions u, v <E C + (D), we say that u is strictly separated (from above) from v if u is 
compactly supported and u{x) < v(x) for every x € supp(u), in this case we write u -< v. 

Definition 2.4. Given u e C + (Q), we say u is a nonnegative viscosity supersolution of (|1.3|) if and 
only if the following conditions are satisfied: 

(i) For every Pq £ Q where u(Pq) > and every function <f> £ C 2,1 (Q) which touches u from below at P$ 
we have at the point Pq 

(Jh >k<j) Aoo^+li^l 2 ifDcfr^O, 

4>t~^k(t>\{D 2 (t>) ifD(f) = 0. 

(ii) If w is a classical moving free-boundary subsolution of (|1.3|) which is strictly separated from u at time 
t\, w(-,ti) -< u(-,ti) and which satisfies w < u on T tlt2 = {(x,t) £ d p Q | t\ ^ t ^ t 2 } , the portion of 
lateral boundary of Q for times between t\ and t 2 (when the space domain is not the whole M. d ), then 
w can not cross u for times in [ii,^]? *-e. w(-,t) ^ u(-,t) for every t £ [ii,i2]- 



See |BV| for a motivation of these definitions. The main idea of comparing with free-boundary solutions 
goes back to [CVj . Finally we can define viscosity solutions. 

Definition 2.5. A function in C + (Q) is a nonnegative viscosity solution of (|1.3|) if it is simultaneously 
a viscosity subsolution and a viscosity supersolution. 



Other ways to define viscosity solutions 

It is possible to define viscosity solutions directly for equation (jl.lj) 

(2.2) p t = Ao. p m 
or, taking w = p m and f3(z) = z~ , for the alternate equation 

(2.3) f3(w) t =A 0o w. 

This last formulation, in particular, has exactly the same kind of degeneracy as the equation (|1.3p we use 
to define viscosity solutions and obtain existence. However, the "Barenblatt" solutions for these equations 
would not have a nonzero normal boundary derivative, and hence they seem to be less natural regarding the 
free boundary propagation. It is interesting to note, however that using the transformations u — -^h[P m ^ 1 
and w = p m on test functions, we can go back and forth from the definition of viscosity solution of one 
formulation to the other. In fact, without having defined it, we will use this "jumping" between formulations 
in the proof of Theorem 17.31 



3 Strictly positive solutions on bounded domains 

To obtain maximal viscosity solution for (|1.3[) . we have to approximate the solution from above by positive 
solutions. Hence, we need to establish comparison and obtain estimates for positive solutions. In what 
follows, Q denotes a cylinder of the form il x (0, T), where fl is an open and bounded domain in K d , and T 
its parabolic boundary, r = (ft x {0}) U (dCl x [0,T]). 

If a function u is a strictly positive classical solution of (|1.3[) in Q, in the sense that there exists c such that 



u(x) ^ c > for every x € Q, 
then u is also a viscosity solution (with obvious adaptations in the above definitions) of the equation 

(3.1) u t = k(3 c (u) Aooii+ \Du\ 2 . 

where /3 C is as defined in the notations at the end of the introduction. Therefore, it is enough to establish 
the results for such functions. 

3.1 Comparison 

We prove two comparison results. The proof of the first comparison result is typical and its idea is similar 
to the proof of Lemma [2~2l however it is weaker than the second. We sketch it here for convenience. 

Theorem 3.1. Assume u and v are viscosity sub- and super solutions of (|3.1[) . respectively, in Q and satisfy 
u,v Js c, 

u(P) < v(P) 

for every P G T . Then u ^ v in Q. If v is a strict supersolution or u is a strict subsolution, then this 
inequality is strict. 

Proof. Step 1. Let us assume the statement of the theorem is false, 

sup (u(P) - v(P)) > 0. 
PeQ 



If we consider instead of v the function v(x, t) = e 7t v(x, h(t)), where h(t) = -, with 7 sufficiently small 

we must still have v > u on d p (tt x [0, T — £(7))) (where £(7) is a small number depending on 7) and now v 
is a strict supersolution. To see this, assuming v is smooth, we compute 



Vi 



jv + e 2 ' yt vt(x, r) ^ 7« + kv A^ v + \Dv\' 



This computation can be carried for the test functions in the definition of viscosity solution, whence our 
claim follows. Therefore we assume v is a strict supersolution. By considering a shorter time interval we can 
also assume 

sup (u(P) - v(P)) =0. 

This supremum must occur at an interior point, Pq = (xo,to) G Q- 
As in the proof of Lemma 12.21 we take 

Wj(x,t,y,s) := u(x,t) -v(y,s) - j\x-y\ - -{t-sf 

and let Pj = (xj , ij , j/j , Sj ) be a point of maximum of Wj inQxQ. Once again, for j large enough, Xj ,yj G 0, 
t/, Sj G (0, T) and P^ — >• (xo, to, xo, to) as j — s- 00. 

Step 2. If Xj = yj, then, with 6j as in Lemma 12.21 (v — Oj)(y, s) has a local minimum at the point (yj, Sj), 
and since v is a strict supersolution we get S < j(tj — Sj). Similarly, u — Q has a local maximum at (xj,tj) 
from where we get the contradiction ^ j(tj — Sj). 

Step 3. Since u(Po) = v(Po), the case Xj 7^ yj can be treated exactly as in the proof of Lemma \2. 2 1 

If either u is a strict subsolution or v is a strict supersolution, then we can obtain the contradiction without 
rescaling the function, and hence, without taking the limit of the rescaled solution and loosing the strict 
inequality. □ 

Corollary 3.2. If u is a subsolution and u ^ M on T, then u ^ M in Q. Similarly, if u is a supersolution 
and u ^ c > on T , then u ^ c in Q. 

Proof. These facts follow from the theorem noting that M + s and c — s are positive solutions of the equation 
for sufficiently small e. □ 



The comparison theorem 13.11 is too weak to prove uniqueness because we assume a strict relation v > u on 
the boundary. However, the transformation we used in its proof to obtain a strict supersolution v above the 
subsolution u will work with an extra Lipschitz condition on the data. 



Theorem 3.3. Assume u and v are viscosity sub- and supersolutions of (|3.1j) in Q, respectively, they satisfy 
u, v ^ c > 0, and 

u(P) ^ v(P) 

for every P G L. Assume also that either v t is bounded below or u t is bounded above on the lateral boundary 
of Q . Then u ^ v in Q. 

Proof. Step 1. The proof is the same as for Theorem 13.11 once we check that for 7 sufficiently small the 
function 

e 7*+7 3 _ 1 

v{x, t) = h'(t)v (x, h(t)) , with h(t) = , 

7 

stays strictly above u on the parabolic boundary of Q. This is clear on the bottom portion; we claim this 
is also the case on the lateral portion of the boundary. We give the proof under the assumption that v t is 
bounded below. The argument when u t is bounded above is completely symmetric. 



Step 2. Let us assume for the moment that the following condition holds, 

(3.2) vt > ~. 

Let us fix x G dQ and write g(t) = v(x,t). Suppose that for some t > 

v(x,t) = h'{t)g(h(t)) < g(t) = v(x,t). 



Noting that 
and 

we must have 



h'{t) - 1 = e 7t+7 - 1 = 7* + 0( 7 2 ) 



t 2 



h{t)-t=^- + 0( 7 2 ), 



g(h(t))-g(t) ,, h'(t)-l /N /N 2o(t) 

• (*) = g ^ ( ;;_f + o(D < -j^g® + o(D = -^i + o(i). 

This is in contradiction with (|3.2|1 for 7 sufficiently small, therefore we must have v(x, t) > v(x, t) as claimed. 

Step 3. Finally, we remove the condition (|3.2[) . Since we are assuming vt bounded below, and v ^ c > 0, 
this condition must surely hold for small f, hence, v ^ u on Q X [0,t] for small i. Since the equation is 
invariant under translations in t, we can now reapply the result to the intervals of the form [kt, (k + l)t], 
k = 1,2,3,.. ., successively, considering u(x, t) = w(x, t — kt) and v(x, t) = v(x, t — kt). D 

4 Special solutions 

We look briefly at concrete examples of solutions with a free boundary which are smooth on their set of 
positivity. All the solutions obtained in this section are classical moving free-boundary solutions in the sense 
of Definition 12.31 These solutions are important as particular examples, and also as models of asymptotic 
behavior of general classes of solutions as t — >• 00 . 

4.1 Separation of variables 

Solutions of the form 

p(x,t)=T(t)F(x) 

are usually of interest. Using the equation, we get the explicit formula for T 

(4.1) r(i) = [C+(m-l)At]~^ T 
and the equation for F 

(4.2) A OOJ F im (x) + AF(a;)=0. 

A special case of separation of variables is when p is independent of t, which corresponds to taking A = 0. 
In this case the equation is nothing more than the infinity Laplace equation. Thus, if v is any nonncgativc 
solution of Aoo v — in fi, then p(x, t) := v~ (x) is a solution of (jl.ip in Q. 

Solving the nonlinear elliptic equation (|4.2[) is nontrivial. This can nonetheless be done explicitly when F 
is radial (we show there is a solution in a non-radial domain in Theorem 17. 3[) . Since the solutions obtained 
this way are known to exist, we postpone their derivation to Appendix A (although their explicit form seems 
to be new) . As observed in the Introduction, the radial solutions are essentially solutions of the 1-d (PME) . 
In fact, all the explicit solutions we obtain are of this kind. 



4.2 Similarity solutions 

If we look for nonncgativc solutions of (|1.1[) of the form 

p(x,t) = t- a f( V ), 
with rj = £~^x and further assume that / is radial and a = /3 = l/(m + 1) we obtain the following solution: 



f(r) 



P{m- 1) 
— i 

2?7l 



By analogy with the porous medium equation, we call this type of solution the Barenblatt functions. The 
solution p takes the form 



(4.3) 



MM) 



Tn 



t^T 



(m^) 2 -\xf 



where j m = [(m — l)/(2m(m + 1))] m_1 . Observe that R denotes the radius of the support of /3r at time 
t = 1. As observed above, these Barenblatt functions are closely related to the Barenblatt functions in 
BV, formula (4.2)], more precisely, they are the (m — l)-th power of our functions in dimension 1, which is 
consistent with the intuition that Aoo is a one dimensional second derivative. 

We will need these solutions for the modified problem, that is, we want to apply the pressure-to-density 
transformation u = (m/(m — l))p m_1 . Therefore, we get the Barenblatt functions for (|1.3[) 



B R (x,t) 



1 



(m + l)t 



Hi- 



lt will, of course, be convenient to consider translations of these functions, in particular in t, to avoid a 
singularity at t = 0: for (x ,t ) S R d x Rq", define 

Sr(x, t; xo, t Q ) := B R (x - x ,t + t ) 

Remark. The choice a = (3 is the only one that makes j md [p(x, t)] d dx constant in time. This is in some 
sense the "correct" behaviour for the infinity-Laplacian, as the "fundamental" solution of u t = A^ u, 
r(x,t) — t~2 cxp(— |x| 2 /4i) satisfies this type of conservation. 

4.3 Traveling waves 

For functions which depend only on one space variable, the infinity-Laplacian is just a regular second deriva- 
tive. Therefore, if we take p of the form p{x, t) = f{rj), r\ = x\ — ct, with c > 0, we have 

(/ m )"+c/' = 0. 

This is the exact same equation obtained, under similar assumptions, for the regular porous medium equation 
(PME), u t = A(u m ), so we get the same solutions, namely, 

p(x, t) = c[a + ct — xi]+. 

The traveling wave solution for (|1 .3[) is therefore given by 



u(x,t) 



m — 1 



c\a + ct — x 



U+ 



For more general planar waves of the form p(x,t) = A(t)U(xi — s(t)), the equation for A and U is again 
the same as for the (PME), and we still have the same solutions. 

It is easy to see that both the Barenblatt functions and the traveling waves satisfy the definition of classical 
free boundary solutions. 
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5 Maximal viscosity solutions in a bounded domain 

In this section we prove existence of the maximal viscosity solution of p.3[) for nonnegativc data as stated 
in Theorem 11.21 The idea is to approximate the data with a strictly positive function. 

Proof of Theorem \1.S\ Step 1. For n e N let g n = g + — and let u n be the unique viscosity solution of (11.31) . 
satisfying u n = g n on T, given by Theorem 11.11 From the comparison results, we have that < u n ^ u l in 
Q for I < n. Hence, (w n )„ e N has a limit as n — »• oo, u. From the estimate of Thcorem ll.il we immediately 
see that u is continuous — if g is Lipschitz, we see that u is also Lipschitz continuous. Note also that because 
u is continuous and («")„ is decreasing, the convergence is locally uniform. 

Step 2. Next we check that u is a viscosity subsolution of (|1.3|) . The argument goes exactly as in |BV( 
Proposition 5.1] for the nonsingular part of the inequality in Definition 12.11 Take a smooth function <fi 
touching u from above at a point Pq = (xo,to) and assume first that D<fi(Po) ^ 0. Since u n converges locally 
uniformly to u, in some neighborhood of Po, Q r = B r (xo) X (to — r 2 , to], for large n the function </>"=(/> + c n 
touches u n from above at a point P n — > Pq as n — > oo. Note that c n is chosen so that the maximum of 
u n — (j> n (at P n ) is exactly and for n sufficiently large D<f) n ^ in Q r . It is clear that c„ — > as n — >• oo, 
and since u n is a positive viscosity solution, 

<#> 0<TAoo</>" + |WT atP„, 
letting n — > oo we obtain the required inequality for cj>, 

If D(j)(Po) = 0, from Lemma [2~2l we can also assume D 2 (f>(Po) = and need only check that <f) t ^ 0. With a 
similarly constructed family (0") ne pj, we have that at the point P n either 

$ < k(f, n Aoo 0" + \Dcp n \ 2 , if D0"(P„) ^ 0, or 
$ < k4> n K{D 2 cj> n ), if D<^. n (P„) = 0. 

In any case, letting n — > oo, we get 

as required. 

Step 3. To prove that u is a viscosity supersolution we need to check conditions (i) and (ii) of Definition ^. 41 
Condition (i) follows just as in step 2 of this proof. Let us check the second condition. Assume v is a classical 
moving free boundary subsolution of (|1.3[) strictly separated from u at t = t\ and on the portion of the lateral 
boundary T tl .t 2 - From our construction, we have u ^ u n in Q, in particular v < u n at t — t\ and on T tl _ t2 - 
Since v is classical at points where it is positive, we obtain from the comparison result for positive solutions 
that v can not touch u n from below. Passing to the limit we get that v can not cross u. 

Step 4. Finally, let us check that u is maximal. Suppose w is another viscosity solution with the same 
data as u which is strictly greater than u at some point. For n sufficiently large, w touches u n from below 
for the first time at some point P n = (x n ,t n ) with t n > 0. Since w = u = g on V and u n = g + — 
on this set, P n must be an interior point. Furthermore, since w(P n ) = u n (P n ) > 0, it is possible to find 
Q n cfix [t n — 5 n , t n ] such that P n S Q n , w > in Q n and w < u n on d p Q n . By the comparison result, we 

get a contradiction. To sec that Q n must exist with the desired properties, note that since w = g = u n 

on r and w(P n ) = u n (P n ) > 0, the open set W n = {x G O | w(x,t n ) > mmw(P n )/2} contains the point P n 
and w = w(P n )/2 on dW n . By continuity, there exists 8 n > such that w > w(P n )/A on W n x [t n — 5 n ,t n ]- 
This set satisfies the above requirements. D 

Maximal solutions are ordered according to their data. 
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Theorem 5.1. If gi and gi are functions in C(T) and g\ ^ gi, then the maximal solutions u\ and ui, 
obtained by the above process with data g\ and gi, respectively, satisfy U\ ^ Ui in Q. In particular, maximal 
solutions are unique. 

Proof. The approximations g" and gj from the previous proof satisfy g" $S g% , therefore, from the strong 
comparison result we have w™ ^ 11%, and in the limit u\ ^ ui. D 

It is worth noting that in Step 3 of the proof of existence of a maximal solution, we in fact prove that a 
maximal solution which is strictly separated at some time from a solution below can not be crossed by this 
solution at later times. 

Theorem 5.2. Let u be a maximal solution and v another solution such that v -< u at t = t\. Then v ^ u 
for all t > t\. 

Here u -< v means u is strictly separated from v as defined in Section [2] The following facts are useful and 
easy to check. 

Theorem 5.3. The Barenblatt functions and the traveling wave solutions of the previous section are maximal 
solutions (in the whole space) of ()1.3[) . 

Proof. We check this for the Barenblatt functions of the form Bn{x,t + to) with to > 0. Take r 3> R and 
consider the approximation given in the proof of Theorem 1 1.2 1 at the beginning of this section, u n with initial 
condition g n = max{l/n, B(-, to)} on the cylinder Q r = B(0, r) x [0, t]. In radial coordinates, the equation is 
simply the regular (PME) in one spacial dimension, therefore, as was proved in jBVj . as n — > 00, the solution 
converges to the Barenblatt function in Q r . A similar argument works for the traveling wave solutions. □ 

It is worth recording the following fact about the independence of the maximal solution on the approximation 
procedure to obtain it. 

Theorem 5.4. Let {h n | n ^ 1} be a family of functions defined on T satisfying h n ^ h n+1 > 0, h n is 
continuous on T and Lipschitz continuous on the vertical portion ofT for each n ^ 1 and lim n h n = g. Let 
v n be the unique solution of (|1.3|) . Then lim„ v n is the maximal viscosity solution from Theorem \1.2i 

Proof. The proof of Theorem 1 1.21 can be reproduced to show that v = lim„ v n is a viscosity solution of (|1.3p . 
It is therefore enough to observe that for any given n ^ 1 there exists m > n such that g m < h n on T (where 
g n is as in the proof of Theorem II. 2|) and use the comparison result to conclude that u m ^ v n . Now take 
the limit in m and in n to see that u ^ v. □ 

Uniqueness of viscosity solutions. If the domain £1 is sufficiently nice and the normal derivative on d£l 
does not vanish, then we can obtain a uniqueness result. We will assume £1 satisfies the following contracting 
property, which is slightly stronger than being star shaped: there exist xq £ f2, A* < 1 and a > such that 

(CP) for every x <G <9f2 and every A G [A*, 1), d(xo + X(x — Xo),dfl) ^ aA. 

Theorem 5.5. Assume dft E C 1 and il satisfies the above contracting property (|CP[) . Assume further that 
uq G M /1 ' 00 (J7) and that the normal derivative of uq on dfl satisfies 



dug 
dv 



sC -m < 0. 

an 



Then the maximal solution of the Dirichlet problem (|7.2[) is the unique viscosity solution of this problem. 
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Proof. Suppose u is the maximal viscosity solution and v is another solution of (|7.2p . After a translation we 
can assume xq = 0. For A E [A*, 1) define u x : ft\ — >• K, 



u x 



(*,*):= A 2+ ^(|,A^), 



where £l\ := Ail = {Ax | x € £1} and 7 will be chosen below. Likewise define Uq(x) = A 2+7 ito(x/A). The 
function u x is a viscosity solution of the equation in fix C fl with initial data Mq . From the assumptions on 
H, and Z?mq we see that choosing A ^ A* sufficiently close to 1 and for x € SI \ fl^ we have 






= (2 + 7)uo(x) — Duq{x) ■ x ^ ma > 0. 

A=l- 



On the other hand, for x <G O^ we can choose 7 independent of A > A such that 
d 



d ^{x) 



(2 + 7)u (x) - Du (x) • x ^ (2 + 7) min u - i-R > 0, 



a=i- fi > 



where L is the Lipschitz constant of uq and i? is the radius of a ball containing £7. Hence, we conclude that 
for A € (A, 1) the function Ug is strictly below uo in tt\. 

We can apply the comparison result in Q\ := Vi\ x [0,T] to conclude that u A < i> in this set. Letting 
A — > 1 _ we obtain u ^ v. Since u is maximal, this means that u = v. □ 

6 Existence and regularity. Proof of Theorem 11.11 

The existence portion of Theorem 1 1.1 1 is obtained, with the necessary adaptations, following the approxima- 
tion procedure introduced in jJKj . Define, for e ^ and 5 > 0, 

C^ 5 u = e A u + , ^"^ (D 2 u ■ Du, Du) , 

\Du\ z + o* 

with /3 C as above, and consider the equation 

( , U t =C^ s f + \Df\ 2 inQ, 

( ' ) \f(P)=g(P) onL 

We estimate the solutions of the approximations and then first let e — > and then S — > 0. 

6.1 Lipschitz estimate in time 

We start with Lipschitz regularity in £ and then prove the regularity in x. 

Theorem 6.1. Suppose g e C 2 (Q), g ^ c > and / = f e - is a smooth solution of (|6.ip . TTien i/iere exists 
K\ > depending only on H-D^lloo, |jD<7||oo, ||<7||oo and ||<?t||oo such that 

\f(x,t)-g(x,0)\^K 1 t 

in Q. If g is only continuous in x and bounded in t then the modulus of continuity of f on CI x [0,4*] (for 
small t*) can be estimated in terms of \\g\\oc and the modulus of continuity of go '■= ff( - ,0) in x. 

Proof. Step 1. Assume g £ C 2 and A > 0. If we define 

v+(x, t) = g Q (x) + A(e At - 1) = g(x, 0) + A(e At - 1), 
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then v + is a supersolution. Indeed, 

v+ -C e - S v+ -\Dv+\ 2 



J t 



= AV< - e A 90 - k ^±^-^- (D 2 g ■ Dg , Dg ) \Dg f 

> e xt (A 2 - p 2 5o||ooA) - ((e + fc||. 9 o||co)||# 2 So!|oo + || DgofJ 
>0 



if we choose A > large enough. In fact, if ||-D 2 .go||oo ^ 1 we can choose A > 1 + y/{l + k\\g\\ oo) + ||-D 2 go||oo 
and if ||-D 2 ffo|| > 1 we can take A > max(2||D 2 g ||, 1 + fc||go||oo + ||A9o||D- Also, for x e dfl and t > 0, 

v + (x,t) = g (x) + X(e xt - 1) >go(x)+t\\g t \\ 00 > g(x,t), 

if A is large, for example A > y 1 1 gt \ | oo • Clearly, v + (x, 0) = g(x, 0), therefore, by the classical comparison 
(see for example |LSUj ) . 

f(x, t) < v+(x, t) = ff0 (x) + A(e A * - 1), 

for every (x,t) £ Q- Similarly, with 

v- = g (x) - A(e At - 1) 

we obtain the symmetric inequality, and hence 

(6.2) \f(x,t)-g (x)\^K t 

for t £ [0, T] where Ao is a constant depending only on the stated norms of g (as it is, it also depends on T, 
but we can take this inequality on a bounded interval [0,9] and then iterate it). 

Step 2. Assume now go is only continuous and let uq be its modulus of continuity. Let us fix a point xq £ il 
and < p < min(dist(xo, dfl), 2\/\\g\\oo)- Let us also define 

g ± (x,t) = g (x )±ujo(p) ± ^f^l* " x ^' 

P 

It is easy to see that g~ ^ g ^ g + on L and thus, again from the comparison principle, f~ ^ / ^ / + , where 
/ ± is the solution of (|6.1[) with initial and boundary condition <jF. Since g^ are in C 2 (]R d x R), we can use 
estimate (|6.2[) to conclude that 

\f ± (x ,t)-g^(x )\^K+t 

where K^ depends on \\g\\oo and p. Therefore, 

\f(x Q ,t) -g (x )\ 

< \f(xo,t) - ^(xo,*)! + l/^so,*) - 5 ± (xo,0)| + 1^(^0,0) - g (x )\ 

< ^l/ + (zo, t) - /-(^o,t)| + A- +i + wo(p) 

< ^\.f + (xo,t)-g + (x o ,0)\ + l\f-(x ,t)-g-(x ,0)\ 

+ -\g + (xo,0) - g-(xo, 0)| + A+t + w (p) 
^2A+i+^o(p). 
With this inequality it is straightforward to conclude the proof. □ 

The full Lipschitz estimate in time now follows easily. 
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Theorem 6.2. If f is a solution if (|6.ip in Q and g £ C 2 (Q), then there exists K 2 > depending only on 
||-D 2 ff||oo, H-Dflilloo, Halloo and \\g t \\oo such that 

\f(x,t)-f(x,s)\^K 2 \t-s\ 

for every x £ Q, t, s £ (0, T). If g is merely continuous, we can estimate the modulus of continuity of f on 
Q in terms of |j<?||oo and the modulus of continuity of g. 

Proof. Taking r > and 

f(x,t):=f(x,t + r), 

using the lemma it is immediate to get 

\f(x,t)-f(x,t)\^K 2 t 
in Ox [0,T — r]. The case when g is only continuous is done in a similar fashion to the previous proof. □ 

6.2 Holder continuity in space 

Theorem 6.3. Let f be the strictly positive solution of (|6.1[) with g £ C 2 {Q) fl Lip(Q), g ^ c > 0. There 
exist a £ (0, 1), and constants K3 > 1 and p 3 , depending only on a, ||g||oo, || A^lloo, ||<7i||oo an d a, such that 
for every e and S sufficiently small and for every Pq = (xo, io) G T and x £ $7 with \x — xq\ ^ pz we have 

\f(x,t ) -g(xo,to)\ ^ K 3 \x-x Q \ a . 

Proof. Step 1. Let us define 

v+(x, t) = g(xo, to) + K*\x - x \ a + A (t - t) , 

where K* Jj 1 and A > arc constants which we will choose in such a way as to make v + a supersolution 
lying above / on the appropriate domain. Let us take x £ fin B P3 (xq) and t £ (0,t ). An easy computation, 
using the fact that \x — Xq\ ^ p 3 ^ 1, yields 

v+ - C e ' s v + - \Dv+\ 2 = -A - K 2 a 2 \x - x \ 2a ~ 2 

,„ ■? I (1 — a)kv + (x,t) ., . \ 

K*a\x - x \ a ~ 2 [ - — K — ; m \]^ j -e(d+a-2)\ 



l + (<J|a:-xo| 1_£ 7#*a) 
K*a I (1 — a)kc 
\x -xo| 2 -« I l + (5/K*aY 



> -^ + — — \2- a I , r , T , — 72 -e(d + «-2) - if*a|«- a; | 



Now note that if J ^ a < 1 then 1 + (5/K*a) ^ 2; furthermore, we have e(d + a — 2) ^ (1 — a)fcc/4 for 

1 — a)kc 
i(d-l) 



e ^ 4(J-i) ' * ne inequality being trivial if d = 1 . We can also choose p% < 1 such that 



ff,a|a;-a;o| < - —!— 

o 

for \x — xq\ ^ P3. Indeed, this is the case if we require that 

aK*p% ^ (1 - a)fcc/8. 
Hence, for 5 and e in the specified range and with this choice for P3, K* and a, we have 

v+ - C s < s v + - \Dv+\ 2 > -A + /<r.apf- 2 (1 ~ a)fcc > 

8 
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provided 

f a-2 ^ 8A 

(1 — a)kc 

Step 2. We want to prove that v + > / on Q* = (12 D P P3 ) x (to — i*,io)i where we take £* := min{l,£o}- 
Let P = (x, i) € d p Q* . Let us first assume P is on the lateral boundary of Q* . If x G 9f2, then, since / = g 
on <9fi and ps < 1, 

/(P) ^ /(P ) + ||P.g||oo|x - x | + ||5t||oo(to - *) 

< <?(P ) + #*|x - xo| Q + A(t -t)= v+(P), 
provided K* Js HDgHoo and A ^ ||gt||oo- If, on the other hand, x € f2ndP P3 (xo), then, using comparison, 

/(P) < II.9II00 < /(P.) + a> 3 + A(t - 1) < «+(P), 

provided K* ^ |MU/P3- 

Step 3. We have to consider the case when P is on the bottom of the cylinder Q* . Let us first assume 
x S ft fl P P3 (xo) and £ = to — 1. In this case, again using comparison, we get 

/(P) < Hfflloo < /(Po) + K*\x - x | Q + A = «+(P), 

as long as A ^ ||ff||oo- Finally, when to < 1: an d hence Q* = (CI (1 B P3 ) x (0,to), we have that / = g on the 
bottom, therefore 

f(P) = f{x, 0) = g(x, 0) < <?(x , t ) + \\Dg\Ux - x | + ||<?t|U*o 
< <?(P ) + #*|x - xo| a + At = v+(P), 

provided, once again, K* ^ ||Z)(7||oo and A ^ ||fft||oo- 
Step 4. To summarize, we have v + ^ / on d p Q* , and hence, by comparison, v + ^ f in Q*, if 

Q (l-a)fcc 

A *"P3 < g 1 

J (1 - a)kc 

A > max { \\g t || «,, IMU}, 

IMI: 



X* ^ max < 1 1 Dg \ 



I oo 1 

Pi 

(1 — a)kc 
S^a and £ < -^_. 

This can be achieved, for example, taking 

,., - . Ilfflloo kc\ . (l kc 1 ., ., i 

A = max{ 5 t oo, j M } , p 3 < mm < . , — f , a = mm <^ -, ^ f and K„ = 5 ooft • 

UIA9II00 16 J L2 16||.g||oo J 

Therefore, we have 

f(x,t ) - g{x ,t ) < w + (x,t ) -fl'(Po) = A'*|x-x | Q . 

Using the barrier v~ := <?(Po) — K*\x — xo| Q + A(t — to) we get the reverse inequality, 

f{x,t ) - g(x ,t ) > -K*\x-x \ a . 

□ 



16 



We can extend the estimate to the interior of the domain. We will use the following notation for convenience. 
For z £ R d , we define Cl z = z + Cl — {x + z \ x £ Cl} and for r > 0, Cl r = {x £ Cl | dist(x, dCl) > r}. Given 
x, y £ M. d we define the closed segment [x,y] :~ {8y + (1 — 6)x \ ^ 8 ^ 1}. The semi-open and open 
segments [x,y), (x,y\ and (x,y) are defined analogously. 



Theorem 6.4. The conclusion of Theorem \6.3\ is valid in the interior of Q, that is, there exists K 4 , depending 
only on ||<?||oo; H-D5II00 and ||gt||oo, such that for e and S sufficiently small and for every x, y £ Cl 

\f(x,t)~f(y,t)\^K 4 \x-y\ a . 

Proof. Step 1. Take a vector z £ B P3 (0) and let V = SI n il z . Let us define f z (x) := f{x — z). From 
Theorem 16.31 we have that \f{x) — f z (x)\ ^ K 3 \z\ a on dV (x £ dV implies that x £ dCl or x — z € dCl). 
Hence, using the comparison principle we have that f z {x) — K 3 \z\ a ^ f(x) ^ fz(x) + K 3 \z\ a for x £ V. This 
means that whenever x, y £ Cl x — y nCl or x, y £ Cl v - x r\Cl, with \x— y\ ^ p 3 , we have \f(x) — f(y)\ ^ K 3 \x — y\ a . 
In particular, the same is true whenever x, y £ Cl\ x _ y \. 

Step 2. When \x — y\ > p 3 , using the comparison principle we obtain the conclusion of the theorem taking 
K± = 2||(7|| 00 //?3 . Let us therefore assume that \x — y\ ^ p 3 and x — y <£ Cl\ x _ y \. Let us first further assume 
that [x, y] C Cl. In this case we can take the two segments [x,w] and [w, y], where w — (x + y)/2 is the 
midpoint of [x, y], let z = y — w and note that w, y £ Cl z fl Cl and x, w £ S7_ z n Cl. Hence, from the first step 
of this proof, we have 

l/(s) " f(v)\ < \m - /HI + l/M - f(y)\ < K 3 (\x - w\ a + \w - y\ a ) 
^2 1 ~ a K 3 \x-y\ a . 

Step 3. If the segment [x,y] is not completely in Cl, then we can certainly find Wi,W2 £ dCl D [x,y] (not 
necessarily different) such that [x, w±) £ Cl and (w 2 , y] £ Cl. In this case we can apply Theorem [63] directly to 
get |/(x)-/(wi)| < K 3 \x- Wl \ a and \f(w 2 )-f(y)\ ^K 3 \w 2 -y\ a . Since \f(w 1 )-f(w 2 )\ < \\Dg\\\ Wl -w 2 \ < 
K 3 \wi — w 2 \ a , we easily get the result with K4 = A 1 ~ a K 3 . This finishes the proof. □ 

6.3 Lipschitz estimate in space 

Observe that, even though the function (x,t) i-> K* \x — Xq\ + X(t — t) is a viscosity supersolution of (|I.3[) 
when xq £ dCl — on the backward cylinder, and for appropriate choices of the constants K* and A — it is the 
e Au term that prevents this function from being a viscosity supersolution of (|6.I[) , as it gives rise to a bad 
term with the wrong sign. It is therefore necessary to let e — > to obtain the Lipschitz estimate. We start 
by obtaining a Lipschitz estimate on the boundary. 

Theorem 6.5. Let g £ Lip(Q), g Js c > ; and suppose f is a (strictly positive) solution of 

, 6S] Ut=£ ' S f+\Df\ 2 inQ, 

\f(P)=g(P) onT. 

There exist constants K§ and p§, depending only on ||g||oo, ||-D<?||oo and ||gt||oo (independent of 6 £ (0, 1)), 
such that for every P$ = (xq, to) £ dCl x (0, T) and x £ CI fl B P5 (xq) we have 

|/0Mo) -g(xo,t )\ < K 5 \x-x \. 

Moreover, if g is only continuous, then the modulus of continuity of f can be estimated in terms of ||.g||oo 
and the modulus of continuity of g. 

Proof. Step 1. Let A'*, L*, 8 and A be positive constants and define 

v+{x, t) = g(P ) +L*\x- x Q \ -K*\x- x \ 2 + X(t -t) + 6. 

n 



We will check v + is a viscosity strict supersolution on Q* (Q* defined as in the proof of Thcorem l6.3[) above 
g for an appropriate choice of constants L„, K* and A. Another easy computation gives 

vt-C°' s v + - \Dv+\ 2 

z J^h \ L ^ 2 KJx-Xo\\ 2 

1 | ( s \ 2 

L ^ \L,-2K,\x-x \ ) 

> -A + kcK t - (L* - l) 2 > 0, 

provided \x ~ x \ ^ 1/2K*, L* ^ 2 and K* > ((i* - l) 2 + A)/(fcc). 

Step 2. Wc now need to choose the constants so that v + > / on the parabolic boundary of Q* . Take 
P = (x, t) be a point in T* = d p Q* . If x £ 90, as before 

/(P) = ff (P) < g(P ) + \\Dg\Ux - Xo \ + |MU(*o - *) 
< .g(Po) + (£* - l)|x - x | + A(*o - *) < v + {P), 

provided L* ^ H-Dfl'lloo + 1 arid \x - x \ < 1/2K*. If x £ ft n dB P5 (x ), then 

/(P) ^ Halloo < .9(Po) + (i* — l)|ar — »o| + A(to - t) < u+ ( p )> 

provided L* ^ H3H00 + 1 and yet again \x — Xq\ ^ 1/2K*. 

Step 3. When P = (x, t) is on the bottom of the cylinder Q* , as before we consider two cases. When to ^ 1, 
t = to — 1 and 

/(P)^!l.g||oo^.g(Po) + A< u + (P) 

as long as A > ||<7t||oo, L* ^ 1 and |x — xo| ^ 1/2K*. On the other hand, if to < 1, t = and hence 

/(p) = 5 (P) < lldlU < ff(Po) + A < w + (P) 

under the exact same conditions as for the previous formula. 
Step 4. Hence we have / ^ v + on T* , and thus by comparison on Q* , as long as we take 

A > Ibtlloo, 
L, > max{2, WDgW^ + 1, UglU + 1} , 

v (i*-l) 2 + A 

A * > IT > 

1 
15 2K, 

Therefore, using once more the comparison principle, we have that for x £ O (~l B pz (xo), 

/(s,to) ^v + {x,t ) ^g(P ) + L*\x-xo\ + 9. 

Since 9 is arbitrary, we have 

fix, to) < g(P ) + L„\x - x \- 

Using instead the barriers 

v~(x,t) = g(Po) — i*|a; — xq\ + K*\x — x | 2 + A(< — i ) — 9 
we obtain the reverse inequality, and as a consequence the Lipschitz estimate. 
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Step 5. Let us finally merely assume that g is continuous and let uj g {o~) be a modulus of continuity at Pq. 
More specifically, let uj g be a continuous, decreasing function in a such that \g(P) — g(Po)\ ^ u g {a) whenever 
max{|x — xq\, \t — to\\ ^ a . Let a G (0,to) and define the smooth functions 

g^[x,t) := g(x o ,0)±u g (a)± ~—\x-Xq\ ± i-to . 

If max{|a; — £Co|, \t — toll ^ cr, then 

g~(P) ^ g(P ) - u g {a) ^ g(P) < <?(P ) +w g (cr) < 5 +(P), 
and if max{|x — a;o|, \t — toll ^ er then 

Therefore, if /* are the solutions of ()6.3j) with initial data g ± , by comparison / _ ^ / ^ / + on Q. Since J 1 * 1 
are smooth we can apply the first part of the theorem to deduce that 

\f ± (x,t )-g ± (P )\^K+\x-xol 

where K§ depends on ||g||oo and a. From these inequalities we get 

\f(x,t )-g(P )\ 

< \f(x,t ) - / ± (^,t )| + \f ± (x,t )-g ± (P )\ + \g±(Pb) - g(P )\ 

< i|/+(x, t ) - f'(x, t )\ + K+\x -x \+ u) g (a) 

< \\f + (*M) g + (Po)\ + \\r{xM g-(Po)\ 

l\g + (Po)-g-(Po)\+K\x-x \+uj g (a) 

^2K+\x-xo\ + ^u g (a). 

This finishes the proof. D 

It is now straightforward to obtain the interior Lipschitz estimate. 

Theorem 6.6. Let g and f be as in Theorem \6.5[ For every x,y £ Q and t G (0,T) 

\f(x,t)-f(y,t)\^K 5 \x-y\, 

where K§ is the constant given in that theorem. If g is only continuous, then the modulus of continuity of 
x i— > f{x,t) can be estimated in terms of ||<?||oo and the modulus of continuity of g in x. 

Proof. Step 1. The proof is very similar to the proof of Theorem 16.41 but in this case it is easy to get the 
optimal Lipschitz constant. Take z G R d such that \z\ ^ p$. Define V — il n £l z and let f z (x) := f(x — z). 
From previous theorem we know that |/(a;) — f z (x)\ ^ ^s|z| on dV. Using comparison, we have that 
f z (x) - K 5 \z\ < f(x) ^ f z (x) + K 5 \z\ in V. Therefore, \f(x) - f(y)\ ^ K 5 \x - y\ if x,y € £l x - y , and in 
particular the same is true if x, y G £l\ x _ y i. £l r — {x G O | dist(a;, dfl) > r}. 

Step 2. Suppose now x — y ^ fiia:— j/l an d let us first assume that the whole segment [x, y] = {z G M. d \ z = 
6y+(l — 6)x, ^ 9 ^ 1} is in f2. Let us assume without loss of generality that p = dist(x, dil) ^ dist(y, 90). 
We can find points i,, ^ i ^ n such that x = Xq, x n — y, Xi G [xj_i,Xj+i] (1 ^ i ^ n — 1), and 
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Pi = \xi — Xi-\\ ^ p (1 ^ i ^ n). Noting that Xi,Xi—i £ fl Xi - Xil we can use the previous step to conclude 
that \f(xi) - f{xi)\ < K 5 \xi - Xi-i\, and hence 

n 

\f(x) - f(y)\ < ]T I/O*) - /(a^*-0l < #s|z - y|. 
t=l 

Step 3. If, on the other hand [x,y] £ ft, then we can find points Xi,X2 € dfl n [je, y] such that [x, xi) C SI 
and [x2,y] \ {£2} C ft. We can further choose w\ £ [x,Xi], with |u>i — xi\ ^ P5 and w 2 € [x2,J/], with 
1^2 -x 2 1 ^ ps- Then we apply step 2 above to obtain \f(x)—f(wi)\ ^ K 5 \x — Wi\, \f(w 2 )- f(y)\ < \w 2 -y\, 
while from the previous theorem, \f(wi) — f(xi)\ ^ K§\wi —Xi\. Putting all these inequalities together gives 
the Lipschitz estimate for this last case. 

The proof of the statement with the modulus of continuity follows as in the proof of Theorem IG.4I □ 



We finally prove Theorem ll.il Existence can be proved piecing out the results in theorems 16.2116. 4l and l6.6l 
and using the standard compactness arguments. Uniqueness follows directly from the strong comparison 
result, Theorem 13.31 



Proof of Theorem \l.l[ Step 1. Assume first g G C 2 (Q)nLip(Q). The comparison principle and theorems 
and 16.41 imply that the family of functions {/ ■ } is uniformly bounded and equicontinuous, therefore, for 
some sequence e^ — > 0, f £k ' S — > f s , which, by a standard argument of viscosity solutions (see the proof of 
Theorem 11.21 above for a similar argument), is a solution of (|6.3[) . 

Using now theorems 16.21 and 16.61 (note that the estimate in Theorem 16.21 is independent of e and 5), we can 
use the same compactness argument to find a sequence f k — > f. The usual arguments for viscosity solutions 
work here to show that viscosity subsolution of (|1.3|) and satisfies condition (i) of Definition 12.41 (we make 
the details explicit in the similar argument for Theorem II. 21 below). The second condition is trivial because 
our functions are strictly positive. 

Step 2. Let us suppose that v is a classical moving free-boundary subsolution of (|1.3|) which is strictly 
separated from / at t\ and on the portion of the lateral boundary of Q with times between t\ and t 2 , r tl]t2 
(see Definition 12. 4j) . Suppose that v crosses / for the first time at an interior point of Qt 1: t 2 = QnR d x [t\, t 2 ), 
-P* = (x*,£*). Since f{P*) > 0, we get a contradiction from the comparison result, Theorem 13. 1[ applied to 
a sufficiently small cylinder around P„, choosing c appropriately. 

Step 3. Since the estimates for the Lipschitz constant (or modulus of continuity) of the approximation / is 
independent of <5, and the convergence is locally uniform, the statements of the theorem concerning Lipschitz 
continuity (respectively modulus of continuity) readily follow. □ 

7 The Dirichlet problem in a bounded domain. Large time be- 
haviour 

If D, is bounded but not a ball, or the data are not radially symmetric, then the exact asymptotic behaviour 
is not the same as in the standard (PME). We treat the general situation here, where ft C M. d is bounded 
and 

(7-1) 9o(x,t)= < 

I uo{x) it t = 0, 

and uq is a Lipschitz continuous function, positive in ft and vanishing on dfl. The Dirichlet problem in 
Q = ft x [0, T] is 

\u solves (|1.3|) in Q and 

^'' \u = g onT = dQx [0,T]UQx {t = Q}. 
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Here is a first result on asymptotic behaviour. 

Proposition 7.1. Let u be the maximal solution of (|7.2[) . Then u decays like 0{t~ l ). More precisely, we 
have 

[Fxix)}™- 1 ^ tu(x,t) < ^(x)}" 1 - 1 

where Fi(x) is the profile given by (|A.2j) for two different radii Ri, i = 1,2. 

Proof. Take xo £ f2 and let B\ = Bh 1 (xq) and B^ — Br 2 (xq) be such that B\ (<= fl C Bi- Let u\ and w,2 be 
the solutions given by (|A.4[) with to — ti and R = R4, i = 1, 2, respectively, and choosing ti < such that 
ui < uo < u 2 at £ = 0. The theorem follows from the comparison result. □ 

Remark. If Q is a ball, then we can successively take the balls B\ and B% sandwiching $7 and obtain the 
exact asymptotic behaviour u(x,t) ~ t^ 1 [F(x)] rn ^ 1 , where F is the profile for the ball f2 given in (IA.2I) . 

In order to get a more precise asymptotic behaviour we need to solve the nonlinear elliptic problem 

(7.3) - Ao. F m = \F 

in Q. We will obtain the result backwards, that is, we will prove that there is an asymptotic profile and then 
show that it solves the elliptic problem (|T.3[) with null boundary conditions. 

In order to prove that in fact there is a profile, we need the following lemma which is adapted from a similar 
result for the (PME) of Benilan-Crandall }BC) (see (Vj Lemma 8.1] for a simple proof for the (PME)). 



Lemma 7.2. The maximal viscosity solution of the Dirichlet problem (|7.2[) , satisfies 
(7.4) u ^-~f 

Proof. Let us for convenience define for A ^ 1 the operator 7a : C(Q) i-> C(Q\) 1 Q\ = X [0, AT], 

T\[u]{x,t) := Xu(x,Xt). 

Consider the approximation u n from the proof of Theorem II .2[ Section and define u" = T\[u n ]. It 
is straightforward to check that u^ solves (|1.3[) with data g% = Xg n (x,Xt). Therefore, for A > 1, using 
comparison, we see that u 1 ^ Jj u n in Q. Hence, from Theorem 15.41 u\ = lim„ u^ is the maximal solution of 
(|7.2[) with data g\ = Xg(x,Xt) and satisfies u\ ^ u in Q. In fact, we must have u\ = 7a [u], since u n — > u 
implies u' l x = Tx[u n ] -> T\[u}. 

Since u\ ^ u and these functions are Lipschitz (from Theorem II. ip . we can compute 



du\(x, t) 
< 



= u(x, t) + tut(x,t), 

A=l + 



dX 

which immediately gives (j7.4[) D 

We are now in position to establish the following: 

Proof of Theorem \1.3\ The above estimates for u and u t suggest the following rescaling: 

i 
(7.5) v(x,t)=V{u){x,t):= [ae (m - 1)r u (i ) e( ro - 1 'j| m_1 , 

where a = (m — l) 2 /m. The bounds t~ 1 u{x, t) ^ C, from Theorem 1 7. II and (|7.4[) in terms of v are 

w(x,r) ^ C and v T ^ 0. 
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It is not hard to check that in fact, at least formally, v solves 
(7.6) v t = A x v m + v. 

The estimates for v now imply that there exists G = Go such that 

lim v(x,t) = G(x). 

r— >oo 

The function Fq from the statement of the theorem is just Fn = a - ™- 1 G = a~ m ~ 1 Ga. D 

Theorem 7.3. The function Gn(x) from the previous proof is a positive viscosity solution of the eigenvalue 
problem (|7.3[) in Q with A = 1 and null boundary condition. 



Proof. Step 1. Let <j> be a smooth function touching G from above at a point xo G Q and assume as usual 
that G — (ft has an absolute maximum at xq. Define the functions T),r] e G C 1 (M) by 

Define also <fi £ (x, r) := <f>(x) + f] e {r) and 

/ r ± i \ m ~ i 

i> £ (x,t) ;=V- 1 {<j> £ ){x,t) = — E (z,log[i^JJ 

Now observe that u(x, r) — e (a;, r) ^ G(x) — 0(x) ^ 0, therefore, there exists c € ^ such that v — {<f) e — c e ) 
has a maximum at a point (x e ,r e ). Since i; T ^ 0, it must be that r e ^ 1/e, as for r < 1/e, £!r (a;,r) < 0. 
Noting that the convergence of v as r — > oo is monotone, the uniform Lipschitz estimates imply that G is 
Lipschitz and the convergence is uniform. Therefore, we must have (x s , r s ) — > (xo, oo) and c e — > as e — > 0. 

Step 2. The transformation V is monotone in the sense that, with r = log^™- 1 ), V(u\){x, r) < V(«2)(£, t) 
whenever u±(x,t) < U2(x,t). Hence, the function tjj £ must touch u from above at the point (x e ,t e ), t e = 
cxp{(m — 1)t £ } and consequently, at this point, 

^ e ,t < ^ e A(D 2 Ve) ifZtyfe = 0. 

This in terms of i/> e translates into 

4>e,r < Aoo 0™ + £ if D^, ^ 0, 

e , T - <fe sC A (i? 2 C) if Dcft e = 

at the point (x e ,T e ). Here we used the fact that when Dw = 0, D 2 (w l ) = lw l ~ 1 D 2 w. Now let e — >• 0. 
Observing that |0 E , T | = |e7y'| ^ 2e and £>0 e = Deft, we get 

< Aoo <ft m + if D0 ^ 0, 

s^ A (L> 2 m ) +0 if D<f> = 0. 

This shows G is a viscosity subsolution of (|T.3[) according to the definition of |CGGj . 

Step 3. To prove G is a viscosity supersolution take now a smooth function eft touching G from below at a 
point xq G Q. Let 

d> e (x, t) := 6{x) - e 5- H -, 

er er 
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where a e is chosen so that <p e touches v from below at some point (x e , t £ ). Let us check that such a choice 
is indeed possible. First observe that for a e = e we get 

(f>e ( x , - J = <p(xo) = G(x ) ^ z> f xo, - 

that is, 4> e crosses or touches v. On the other hand, letting t* be such that t > t* implies v(x, r) > G(x) — e 
(recall that the convergence as r — > oo is uniform), if o € < — (t*) 2 £L v /2, where L v is the Lipschitz constant 
for D, then <f> £ < v. Indeed, for r ^ t* we have 4> e (x,r) < G(x) — e < v(x,t) and for r < t* we compute 
6iT (x,r) = —j — -% > L^, which means we must have v(x,t) > <P(x,t). Therefore, there must exist 
a e <G (— (T*) 2 eL v /2,e) such that </> £ touches v from below at a point (x 6 ,t £ ). Note also that the above 
implies r e > 1/e, hence we have lim e _ > o(x e ,T e ) = (xo,+oo). The rest of the argument goes as in step 2 
above. □ 

8 The Cauchy problem 

In this section we consider briefly the Cauchy problem 

(gl) Ut^kuA^u+lDul 2 , inR d x[0,T], 

iu(x,0) =u (x), for xe R d . 

When u : K d — > K is uniformly continuous, positive and has compact support we can readily apply the 
foregoing theory in a sufficiently large cylinder to obtain existence and uniqueness of a maximal solution, 
which necessarily has compact support as well. On the other hand, for positive data u ^ 0, we can 
obtain the solution by considering first strictly positive data and approximating the problem posing it in 
il r = B 2r (0)x [0,T], 

J ut = kj3 c {u) AooU+ \Du\ 2 , in Q r , 

1 u(x, t) = Uq(x), if (x,t) € d p £l r . 

In principle we could take Uq(x) = Uo(x), but for technical reasons we will have it depend on r. In any 
case, we choose it so that {u^} is uniformly equicontinuous on compact sets, Uq Js c > and its modulus of 
continuity is bounded by the modulus of uq. For each r > 0, let 

if |x| < r, 
«5(a:) - <( M if |z| = 2r, 




-M 



} ifN 



if Ixl = 2r - Xr. 



Theorem 8.1. Given uq : R d — > R uniformly continuous, bounded and such that Uq(x) ^ c > 0, i/iere exists 
a maximal bounded viscosity solution of (|8. 1 [) . TTie modulus of continuity of u can be estimated in terms of 
the modulus of continuity of uq and ||uo||oo- 

Proof. Step 1. Let u r be the unique solution of (|8.2[) given by Theorems 11.11 and 13.31 Fix R > and 
consider the family {u r } r> R. By the estimates of Theorem 11.11 this family is uniformly bounded and 
equicontinuous on Qr. Therefore, using a diagonal argument, we can extract a sequence which converges 
locally uniformly on R d x [0,T] to a continuous function u whose modulus of continuity depends on the 
modulus of continuity of uq and ||uo||oo (and is Lipschitz continuous whenever uo has this property). Note 
also that inf K d u = limhnv inf K d u r ^ c > 0. The proof that u is a viscosity solution goes as in the proof of 
Theorem 11.21 above. 



Step 2. To show that the solution obtained in Step 1 is maximal, let v be another solution. We prove in the 
next lemma that v ^ M. Since u r ^ M on the lateral boundary and Uq ^ uq = v on {t = 0}, by comparison 
we immediately obtain the result. □ 
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In the proof of maximality above, we have used the following version of the technical Lemma 5.2 from |BVj . 
Lemma 8.2. If u is a bounded viscosity solution of (|8.ip such that < c ^ uq ^ M , then u $J M everywhere. 

Proof. We consider the same function as in |BV[ Lemma 5.2] 

V(x,t)=N + b^- t +Xt 

and easily check that V is a classical strict supersolution if 

1 > 2(fc + 2)6 and AT(1 - 2kb) > 2kNb. 

These inequalities are satisfied if we select b small enough and A not too small. In particular, this is true 
for b = e, N = M + e and A = ^^-e = O(e) with e small enough. This function is strictly greater than u 
at t = and, because we are assuming u is bounded, it is also strictly above u for all large \x\ uniformly in 
[0,T]. If V is not strictly larger than u everywhere, there exists Pi = (oJi,ti) G M. d x [0, T] where V touches 
u from above. At this point, from the definition of viscosity solution, we have 

v t (Pi.) < fcy(Pi) a^ v(Pi.) + 1 wop)! 2 . 

This contradicts the fact that V is a strict supersolution in a small parabolic neighbourhood of Pi, hence 
we conclude that u ^ V on M. d x [0, T]. Letting £->0we conclude that u ^ M everywhere. □ 

Using the same ideas as above, we can let c — > to obtain maximal solution for nonnegative data. 

Theorem 8.3. Given uq : R d — > R uniformly continuous, bounded and such that uq(x) ^ 0, there exists a 
maximal viscosity solution of (|8.1[) . The modulus of continuity of u can be estimated in terms of the modulus 
of continuity of uq and ||uo||oo- 

8.1 Large time behaviour 

Considering data Uq with compact support, using the Barcnblatt functions and the traveling waves solutions, 
it is interesting that the ideas of the classical (PME) can be used to obtain properties for the support of 
super- and subsolutions and in particular to obtain the asymptotic behaviour as t — > oo. 

Proposition 8.4. We have the following properties: 

(a) The support of any viscosity supersolution is nondecreasing in time and penetrates the whole space as 
t — > oo. 

(b) The support of any viscosity solution expands in a continuous way. 

In particular, any solution whose initial condition has compact support, has support in a set expanding like 
0(i 7;r + T ) as t ->• oo. 

The proofs of (a) and (b) are identical to the proofs of |BV( Propositions 6.1 and 6.2] and the last statement 
follows from sandwiching the support of the solution between to Barcnblatt functions. Wc can now prove 
our final theorem. 



Proof of Theorem \l.J\ Step 1. The first step is to sandwich p{x,t\) for some t\ > between two 1-d 
Barenblatt solutions with 1-d masses < Mi < M^ (note that Mi = 2 J pi(r,0)dr), and this has been 
explained above. Such bounds will preserved for all times f^t\. 
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Step 2. Next, we copy from the study of asymptotic behaviour of the 1-d PME to define a family of rescaled 
solutions p\ for all A ^ 1 as follows 

(8.3) p x (x, t) = X^p(X^x, Xt). 

It is easy to see that the p\ are maximal viscosity solutions of the Cauchy Problem (with rescaled initial data) . 
The bounds with Barcnblatt solutions imply that on any compact time interval [t±, £2] with < t\ < £2 < 00 
the family p\ is continuous, uniformly bounded, and supported in a uniform ball Bn r (0). 

Step 3. Now we use Aleksandrov's principle, as explained for instance in [CVWj . to show that for a solution 
p{x,t) with initial data po(x) ^ supported in the ball -Br(O) we have for all t ^ and all r > R 

inf p(x.t) = max p(x,t) 

\x\=r \x\=r+2R 

If this is applied to the rescaled solutions, we get for all |x| ^ R\ = RX~ 1 '( m+1 > 

inf p\(x,t) = max p\(x,t) 

\x\=r \x\=r+2R x 

Step 4. We now fix t = 1, A very large, so that R\ ^ e is very small, and define 

pi(r) = inf px {x, 1), p 2 {r) = maxp A (i, 1), 

|a;|— r |x|— r 

We easily verify that ^(r), f>2{f) are nonnegative and radially symmetric functions, both supported in the 
same ball Bft_ t (0), they are nonincreasing as functions of r for r ^ e, and we also have 

P2(r) > pi(r) ^ /9i(r + e) 

for all r ^ e. It is then easy to verify that the 1-d mass of P2M — Pi( r ) is less than Ce. 

Step 5. If pi(r,t) and p2(r t t) are the corresponding solutions of the 1-d PME with initial data at t = 1 
given by /5i(r) and P2(^), respectively, we have for all t ^ 1 

pi(r,t) 5^ /9a(K)*) < /5a (>",*) 

In view of the convergence result for the 1-d PME, cf. [Vj Theorem 18.1], the result follows. □ 

A Radial solutions with separation of variables 

Let us look more closely at solutions of the form 

p(x,t) =T(t)F(x) 
with T and F as in (|4.1[) and (|4.2[) . We further assume that F is radial, F(x) = F(r), with r = \x\. 

A.l Radial solutions with A > 

Let g(r) = A"^F m (r). Then from (jUJ) wc get 

5"+.9 p = 0, 

This equation can be integrated and yields 

9(r) „ ( 1 1,, 1 9 p+1 (r) \ 

—j= 2F1 — — -, - 1 + — — -, ) = ±kr + C, 

sfa \p + 1 2 p + 1 a 
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where 2F1 is the hypcrgcomctric function and C and a are arbitrary constants, a > 0, and k = ,/rrr. Let 



us for convenience define H p (z) := -4- 2F1 I — W, 5, 1 H — tj, ^5— ) , which satisfies H'[z) 

1 
that H p is left-continuous at z = a^ 1 . In fact, we can compute 



Va-z p+1 



Note 



H p (0) = and H p (a^+r) = a^+r 



* > r i + iir 



r(i + ^- 

1 2 ^ p+i 



=: A, 



We have solved the equation for g implicitly, H p (g(r)) = ±kr + G, where k = <^/2/(p + 1). Since we want 
g Js 0, we need to define -ff p on the interval I p = [0,ap+ 1 ]. On this interval, H' > 0, therefore H p has an 
inverse, G p : J p — > I p , where J p = [0, A p ], and G p > 0. Then we have 

(A.I) g(r) = G p (±kr + C). 

Depending on whether we choose the plus or minus sign on this identity, we get two different types of 
solution. 



A. 1.1 Solutions denned on a ball 



Taking the minus sign in (|A.1[) we get g(r) = G p (k(R — r)), which is defined and smooth for r £ [0, R] with 
kR = A p , that is 



R = 



^r(i + ^ T 



r 



\2 ^ p+1 



P+l 



a 2 (m + l) 



The function g is radially decreasing, vanishes at r = R and has a maximum at r = 0, 5(0) = af+ 1 = MR™- 1 , 
where M is a constant depending on m. 

The solution p{x,t) = T(t)F(x) is given by (|3~Tj| and 



(A.2) 
that is, 

(A.3) 



F(x) = A~ 



G, 



2m 

TO + I 



(fl-W) 



p(x,t) 



1 



[( m -l)(t-t )]s=T 

For the transformed equation (|1.3|) we get the solution 



G,, 



2m 
m + 1 



(R-\x\ 



(A.4) 



u(x,£) 



G , , 



2m 



j(«-W) 



(m-l)a(t-to) 

Note that u and p have infinite gradient at the boundary of the domain {|x| = R}. 



A. 1.2 Solutions defined on an annulus 



If we take the plus sign in (|A.1[) instead, we get the solutions defined on an annulus {i?i < r < R2}, where 
R-2 = Ap/k. These solutions are 



(A.5) 



p(x,t) 



[( m -l)(t_io)]^ 



G, 



2m 
m+ 1 



(N--R0 
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for (|1.1[) . and 

(A.6) 



u(x, t) 



{m-l) 2 (t-t Q ) 




m + 1 



(M-iZi) 



for (|1.3[) . These functions vanish on the inner boundary of the annulus, are radially increasing and have 
maxima on the outer boundary equal to the maxima which the solutions of (|A.3|) and (|A.4[) attain at x = 0. 
If — A p /k < R\ ^ 0, then the same formulas define solutions on the annulus {0 < r < R2}, which are 
continuous but not differentiablc at x = 0. 



A. 2 Radial solutions with A < 

For negative A, we let £ = |A| and still assuming that F is radial, we let g(r) = £~ rn - 1 F m (r). Then g solves 

g" - g p = 0. 

Integrating once as in the previous case, we now get 

9' =±K rxi 

y/a + gP +1 Vp+1 

We now need to distinguish three cases: a > 0, a = and a < 0. 

A. 2.1 Case a > 

We compute T p (g(r)) = ±J^r + C, where I p {z) = ^ 2F1 [j^, \\ 1 + ^j, -^-J ■ This function is 

1-to-l from [0, 00) onto itself with I p (0) = 0, I p '(0) = -4=, / p (z) ~ j^— z - ^ - asz^ +00. Let J p denote the 
inverse of I p on this interval. The function J p is also a bijection of [0, +00) onto itself, J p (0) = 0, J P '(0) = \/a, 
J P (z) grows like z 1 -? as z — > +00 and J p (z) grows like z 1 -? — z™- 1 . Since g(r) = J p I ±« /— rj?* + ( 
have two possible types of solution depending on which sign we take. If we choose the plus sign, 



F{r) = £—1 



and 



p(x,t) 



u(x, t) 




[(m.-l)(t Q -t)]^ 



m+l 



(\x\-R) 



m 



(m-l) 2 (£ -i) 




m. + 1 



(\x\-R) 



For R > 0, these functions are defined on {\x\ ^ R} x (— 00, io) with u, p = +00 the top portion of its 
boundary, t = to, u, p — on the lateral boundary, \x\ = R, and again Du,Dp = +00 on the lateral 
boundary. As |x| — > 00, p grows like |x| m - 1 and u grows like |x| 2 . For R = the situation is identical. The 
function is defined in all space but Dp has a singularity at x — 0. Finally, if R < 0, then p is defined for all 
x but the gradient, though bounded near the origin, is discontinuous at x = 0. 

If we take a minus sign in the implicit identity for g, using the same inverse for 7 p , we get functions defined 
on a ball with a finite but discontinuous gradient at the origin. 
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A. 2. 2 Case a = 

The integral for g is jz—g^ E (r) = iv^rr 7 * + ^' therefore 

1-p 



In this case wc have 



(A.7) 



g(r) = 

F(x) = ^ 
p(x,i) 



a/2(p+1) 



m — 1 



C 



=(M-A) 



N /2m(m+ 1) 

(777 -l) 2 \^ r (|x|-i?) 2 
2m(m + l)/ [ (i -<) 



and 



i(x,t) 



i (M-i?) 2 



2(m+l) (*o-t) 

For R > these functions can be defined either on the interior or the exterior of the ball {|a;| ^ R}. For 
R < they are smooth in the whole space minus the origin and for R = the function u is a smooth solution 
of (TO}. 

A.2.3 Case a < 

In this case the integral of g is K p (g) — ±-i/— ttT" + C, where A p : [|o|5Tt,+oo) — > [0, +oo) satisfies 



if '( z ) = -i 

A'p = and K p '(z) 



and A" p (ap+ 1 ) = 0. This function grows at infinity like z 2 and at z = \a\p+ 1 wc have 



/ 



(p+l)(«-|o|P+i) 



If we denote by L p : [0, 



|a| p+1 , +00) the inverse of A p , we 



have L p (0) — \a\ P+ 1 , L p '(z) ~ -y/(p + l)z as z -> 0, while asz-> +00, L p (z) ~ 2 •■-» and L p '(z) ~ z'-p . Wc 
then have 



(A.8 
and 



F(ar) = Z" 



L„ ± 



2m 

m+1 



C 



p(x,t) 



[(m-l)(t -t)] — 



± 



2 m 

771+1 



W + c 



Since L p (z) is like a power z 2 at z = 0, the conclusions regarding p and 77 arc similar to the case a = 0, we 
have functions defined on the interior and on the exterior of {\x\ = r}. 

Comments and open problems 

We do not analyze here the fine behavior nor regularity for the free boundary. This might in fact be a key 
component for a full uniqueness result as it happens with the (PME) (see |CVj ). which is possibly the main 
open problem related to this equation. 

Another idea we do not explore is how to develop a numerical scheme, both for the equation and the free 
boundary. 
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We are planning to pursue the ideas of interpolated models mentioned in the introduction separately. 
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